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V. DESIGN OF THE ANALOG ASPS THROUGH DECOUPLING 
AND POLE PLACEMENT 


5.1 Introduction 

In Chapter II the analog controllers of the ASPS are designed for the control 
of the X, <j)^ and dynamics. The eight eigenvalues of the system are assigned 
so that the overall system has an equivalent bandwidth of 2 Hz, We shall see 
later that this Is not an accurate description of the bandwidth requirements of 
the system. 

The analog controllers consist of eight state-feedback gains with no coupling 
betv/een x, <J)| and ^ 2 ' In other words, the Input U| of the x component Is affected 

f 

only by feedbacks from x and x, the Input O 2 of the component Is affected only 
by feedbacks from /4i| , <j)| and (fij , Similarly, the Input u^ of the <^2 component 
is realized by feedback from the states ((>2 and 1 P 2 . Therefore, the three 

independent controllers are essentially PID controllers. 

The values of the eight feedback gains are determined by using the Brown's 
method for pole-placement. The sampled-data version of the x, ({)j and <f >2 dynamics 
of the ASPS was obtained by inserting sample-and-hol d devices in the three input 
channels. It was shown in Chapter III that using the feedback gains designed 
for the analog system, the sampled-data system is stable for sampling periods 
less than or equal to 0.0075 seconds. However, this sampling period Is still 
considered to be too small to be economical and practical for the ASPS. 

5.2 Design of the Analog ASPS Through Decoupling And Pole Placement 

As it turns out the desired bc^ndwldths requirement of the ASPS is as 
fo 1 1 ows : 

Xj dynamics 0.04 Hz 

(j)j dynamics 


10 Hz 
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dynamics 1 Hz 


Since the <|)j dynamics Is 250 times faster than the x dynamics, It Is 
virtually Impossible to find an equivalent bandwidth of the overall system and 
establish a general eigenvalue requirement for the system. Therefore, It Is 
necessary to decouple the x, i|i| and 1 ^ 2 ^ dynamics through state feedback and 
simultaneously place the poles to realize the desired bandwidths for all three 
system components. The problem Is stated ast 

Given the system ORIGINAL PAGE IS 

x(t) = Ax(t) + Bi^(t) OP POOR QUALITY (5-1) 


where 


x(t) = 


u(t) 



(5-2) 


(5-3) 


A is the 8x8 coefficient matrix and B is the 8x3 input matrix. The state feedback 
is defined by 


u(t) = -Gx(t) 


(5-4) 


where G Is the 3x8 feedback matrix. The elements of G are to be selected so that 
the coefficient matrix of the closed-loop system, A - 8G, has the following form: 
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A - BG 


■A. 


A, 


ORIGINAL S 

OF POOR quality 


3J 


(5-5) 


Aj Is a 2x2 matrix for x 

A^ is a 3x3 matrix for <l>j 

A^ is a 3x3 matrix for ‘l> 2 . 

The eigenvalues of A^ ^ A^ and A^ are so selected that the bandwidths requirements 
are satisfied. 

Referring to fig. 2-1 the A matrix of the ASPS has the following form; 


’21 


51 


81 


0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


0 

®2^ 

1 

0 

0 

0 

""84 


0 

0 

0 

I 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


27 


57 


87 


0 

0 

0 

0 

0 

0 

I 

0 


(5-6) 


The 0 matrix is written as 


B s 


21 


51 


81 


22 


23 


52 *"53 


82 


‘’83 


(5-7) 
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where the a|j's and hjj's In the A and B matrices are given as 

a2, “ -0.61207568 
a^ij = -) ,k 8 hk 272 X lo"^ 

= I .A844272 X lO"^ H 
a^i » -0.31202659 
^54 “ -7.589096 x lo’^ 

= 7.589096 xio’^ » 
agj =« 0.31202659 
^84 “ 7.589096 X 10“^ 

^87 ** -7.6884996 X 10 V 

b2, = 0.58237458 
622 = 0.29688544 
» - 0.29688544 
bjl = 0.29688544 
bj2 » 0.1 51 781 92 
bj2 = -0.1 51 781 92 » -bj2 
bgj = - 0.29688544 
bg2 =* -0.1 51 781 92 
bg 2 » 0.15376999 

It is simple to show that 8 has full rank (rank = 3) « 

Let the feedback matrix be represented by 



"^ll 

912 

^13 

914 

^15 

^16 

^17 

9i8 

G = 

921 

^22 

®23 

924 

^25 

^26 

^27 

928 


331 

932 

®33 

934 

935 

^36 

^37 

^38 


Then BG becomes 
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8G 


0 

(BG) 

0 

0 

(BG) 

0 

0 

(BG) 


21 


51 


81 


where 


0 

(BG) 

0 

0 

(BG) 

0 

0 

(BG) 


22 


52 


82 


(BG) 


0 

(BG) 

0 

0 

(BG) 

0 

0 

(BG) 


23 


53 


83 


0 

(BG)2^ 

0 

0 

(bg) 

0 
0 

(BG) 


5^ 


8^} 


IJ 


i, 


0 

(BG) 

0 

0 

(BG) 

0 

0 

(BG) 


25 


55 


85 


0 

(BG) 

0 

0 

(BG) 

0 

0 

(BG) 


26 


56 


86 


0 

(BG) 

0 

0 

(BG) 

0 

0 

(BG) 


27 


57 


87 


0 

(BG) 

0 

0 

(BG) 

0 

0 

(BG) 


28 


58 


88 




(5-9) 


In order for A - BG to be a diagonal matrix, the following relations must 


nold; 


(BG)23 “ ^23 « 0 


24 ®24 


"25 “ ° 


(BG) 

(BG)25 

^^*^^27 “ "27 

^®^^^28 “ 328 = 0 

^^*^^51 " "51 
(86)52 = 352 “ ° 


(BG) 

(BG) 

(BG) 

(BG) o, = a 


56 " "56 

57 “ "57 

58 " 358 


= 0 


= 0 


(BG) 


81 

82 


=5 a 


81 

82 


= 0 


(BG)8'3 “ "83 “ 6 

^86)34 " "84 
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(BO) 85 * “85 ■ ° 

These equations represent 1$ equations with 24 variables In the feedback gains 
9ij (l"l»2j3l j“ I > 2 1 » . . j 8) « 

For the present design, the x dynamics Is specified to have a bandwidth of 
0,04 Hz. Since the x dynamics are represented by a second-order system, the 
bandwidth of the system Is given by 

BW » 0)„( 1 - ^ - 4?^ + 2 (5-11) 

n ' 


If we choose the damping ratio ? to be 0.707, then 

BW = m (5-12) 

n 

where Is the natural undamped frequency. 

Therefore, for the x dynamics, the eigenvalues are selected to be at 


where 


Pi » ?2 “ 1 -)®l 

^ , Q.04 X 


(5-13) 

(5-14) 


The <j>^ and dynamics are of the third order. However, we can use second- 
order approximations by placing the real root far to the left on the real axis 
In the s-plane. For this case we let the real root of the tj)^ dynamics be placed 
at 

Pj “ 82 = -200 (5-15) 


Then, the complex roots are: 


Pl ^> P5 = ^2 i -'®2 


where 


^2 


22LJLi^« -44.436 
/2 


Similarly, for the dynamics, we let 


(5-16) 

( 5 - 17 ) 


; P6 = ^^3 “ 


(5-18) 
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Then, 


where 


P7» Pa " ^3 t J^3 

27T X 1 I. I f, / 
®3 “ ” ”72 " “4.«i43o 


(5-19) 

( 5 - 20 ) 


The characteristic equations of the three decoupled subsystems are as 
f 0 1 1 ows ; 


X component: (s - 3| - Jaj)(s + 3| + Ja^) 


2 2 
«» s -*2a,s + 2aT ® 0 


( 5 - 21 ) 


components (s - b2) (s - -ja2)(s - Oj + ja^) 

» (s - b 2 ) (s^ - 2 a 2 S + 23^) 

. , = s^ - (23^ + b^)5^ + (2a^ + 2a2b2)s - 2a2b2 = 0 (5-22) 

(|)2 componwnt: (s - b^) (s “ 23 ” Ja 3 )(s “33 + Ja 3 ) 

a (s - bj) (s^ - 23^5 + 2a^) 

a s^ - {2a^ + b2)s^ + (2a^ + 2a2bj)s - 2a2b^ = 0 (5-23) 


If there exists a set of feedback gains such that the decoupling of the 
subsystems x, ‘I’j and can be accomplished, then the decoupled closed- loop 
system will have the Independent characteristic equations of Eqs. (5-21), (5-22) 
and (5-23). 

Let the matrix A of Eq. (5-6) be written as 

(5-2i») 

where Ajj denotes a 2x2 matrix, A 22 and A ^3 are 3x3 matrices. The submatrices 
denoted by X contain elements which are unimportant for the Immediate development. 
In view of Eqs. (5“5) and (5-9), the decoupled x dynamics with state feedback 




A = 


^ ! ^22 I 

“ r - - ' 


1 


X 

A. 


I "33 


can be expressed as 
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•X 

1 


1 

X 

1 

■ Aj 

L’-aJ 


— 1 
C4 
X 


PAGE IS 
OP POOR QOALITy 


where 


A, . A|| - 

“ 0 

0 

<“=>22. 

m 

• 

" 0 1 " 

.^21 - (“'=>21 ■(“'=>22. 

The characteristic equation of A| 

Is 



. 1 

s 


-1 

1 


si - A, 


^21 " ^ 


*s + (80)22^ *■ ^21 ^^*^^21 ** ^ 

Similarly, the decoupled (j>| dynamics with state feedback is described 


X3 


t 

rr\ 

X 


" ^2 


_^ 5 _ 


_’^ 5 _ 


where 


^2 “ ^22 


0 

0 


0 

0 


0 

0 


(BG)j 2 (BG)g^ 


1 

0 


0 

0 

(BG) 


55 


-.{BG)53 a^^ - {BG)^^ 


0 

1 

-(BG) 


55 


The characteristic equation A2 is 


si -A, I = 


s 

0 


( 88)53 “^ 5 ^ ^ ^^*^^54 


0 

-1 

s + (BG) 


55 


(5-25) 

(5-26) 

(5-27) 

by 

(5-28) 

(5-29) 
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" 3^ + (BG) 5 gS^ + (('^ 0 ) 5 /, - + (BG)gj - 0 

For the (|>2 component, the decoupled closed-loop stete equations are 


t 


— — 

^6 


^6 

X 7 

* 

- A 3 

X 7 

_^8_ 


CO 

X 

1 


where 


^3 “ ^33 


0 

0 

(BGI 


86 


0 

0 

(BG) 


87 


0 

0 

(BG) 


0 
0 

^^^*^^86 ^87 

The characteristic equation of Is 


0 

- (BG) 


87 


88 

0 

I 

(BG) 


88 


sl - A, 


s 

0 


0 

■I 


-=87 * ‘““>87 = <““>88 


“ + (88)335“ + ((88)3, - a3,)s + (bg)„3 • 0 


87 87-' 


86 


(5“30) 


(5-31) 


( 5 - 32 ) 


(5-33) 


For pole placement, In order to meet the bandwidth requirements, the corresponding 
coefficients of the characteristic equations in Eqs. (5-27), (5-30) and (5“33) 
must match those of Eqs. (5-21), (5-22) and (5-23), respectively. 

Thus, for pole-placement, 

^ - 2a, s + 2a, ^ = s^ + (BG) 22 S + (BG)j,, - a^, (5“3^«) 

s^ - (2a^ + ,b2)s^ + (2a^ + 2a2b2)s - 2a^b2 

== s^ + (BG)g5S^ + [(BG)^^ - a^Js + (BG)^^ (5-35) 


/ 



( 5 - 36 ) 


- 5 ^ + (BG)g 8 S^ + ((BG)gj, - ag 7 )s + (BG)gg 

Equating the like coefficients In the above three equations, we have 

( 66)22 “ *> 21^)2 *^ 22^22 ’’ 23^32 " ■’^®1 

(BG)2| - b 2 j 9 j, + b 2232 | + 62393, - 2a^ + a2, 

(BG)55 " V 25 * 953^35 “ -^282 + b2) 

" *^51^14 ‘^52^2^1 *'539311 ■ 232 ^ ^^ 2^2 + ^54 

(=0)53 - b5,g,3 + bj^g^j + bjjgjj - -aa^b^ 

" *'8|9|8 ■*’ *'82928 *'83938 " "^283 + 63 ) 

(0G)87 - bg,g,7 + bg2927 + bg^g^^ « 2a^ + 2a3b^ + ag^ 

“ *'8j9j^ + >582928 + bg^g^^ « -2a^b3 

Notice that the 16 constraint equations on the feedback gains In Eq. ( 5 -IO) 
i»re conditions on decoupling of the three subsystems, whereas the 8 constraint 
equations in Eq, (5-37) are the conditions for pole placement. The 24 constraint 
equations contain 24 unknowns In the elements of the feedback matrix G, g ,1=1, 
2,3, j-1 ,2,.. .,8. 

The 24 constraint equations in Eqs. (5-l0) and (5-37) can be written as 


WG « C 


( 5 - 38 ) 


where 



*'21 

‘'22 

''23 

w = 

*'51 

'52 

*’53 


^^81 

*'82 

*'83 _ 

G Is the feedback matrix given 

In Eq. ( 5 - 8 ), 

and C 

are composed of a, , a. 

, b, , 

62, 63, and 

a|j, 1 


The feedback gain matrix G can be solved form Eq. (5-38) if W is nonsinguiar. 


Then , 


G - W’^C 


(5-39) 


For W to be nonsinjiular, It Is necessary and sufficient for B to have full rank. 

For the system parameters given for the ASPS, the feedback matrix Is solved 
from Eq. (5-39), and the results are given as follows: 


g^ I « 36.I*A9 
gj2 “ 212.^ 

» - 9 . 2569812^*0 X 10 ® 
g,/, » -2.5^7162662i* x 10^ 

9|5 =■ -3.3892568 X 10^ 

916 " 0 

917 = 0 
9,8 “ 0 

92, = - 73.35 

922 “ “^'5.^<5Ai* 

922 “ 2.2126125552 x 10® 
a 6.0882526771 x 10^ 

925 “ 8-'012732 x 10® 

926 “ 3.96750301* x 10® 

927 “ ' -091703152 X 10® 

923 = 1 .1*5266^ X 10^ 



932 = 3^96750304 X 10® 
92f, * 1 .091703152 X 10^ 
= 1 .452664 X 10® 

926 " 3.96750304 X 10® 
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9^7 - 1 .091703)02 X 10^ 

g^Q = 1 .^5266^1 X 10^ 

Substituting G Into Eq. (5-5), the elements of the closed-loop A - BG matrix 
are obtained as follows: 



ORIGINAL PAGE IS 
OF POOR QUALITY 
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VI DESIGN OF THE DIGITAL ASPS THROUGH DECOUPLING 
AND POLE PLACEMENT 


In this chapter the digital ASPS Is to be designed through decoupling and 
pole placement. In Chapter IV the digital ASPS was realized by placing sample- 
and-hold units at the three Inputs to the x, and components. The feedback 
gains were the same as those of the analog ASPS. 

The dynamics of the system are still described by 

^(t) « A^(t) + B^(t) (6-1) 

where A, B, x(t) are as defined In Chapter V. The control vector £(t) Is given 

u{t) = -Gx(kT) kT < t < (k+I)T (6-2) 

The solution of Eq, (6-1) over the interval kT £ t <_ (k+l)T Is 


t 

£(t) = (|)(t - kT)x(kT) + / (j)(t - T)B£(T)dt 

kT 

t 


where 


f(|>(t - kT) - / tf)(t - x)dTBG]x(kT) 
kT 


At 


«j>(t) = e‘ 

Substituting t = (kH*l)T Into Eq. (6-3) » we get 


Let 


Then Eq. 


T 

x((k+l)T] = ((|){T) - / <J)(T - T)dxBG)x(kT) 

0 

T 

0(T) = / t|)(T - T)dxB 

0 

(6-5) i 5 s impl if led to 

x((k+l)T) = (c|)(T) - 0(T)G)x(kT) 


(6-3) 

(6-i|) 


(6-5) 

( 6 - 6 ) 

(6-7) 


Equation ( 6 - 7 ) now represents the difference state equations of the closed-loop 
* 

state-feedback ASPS. 


The design objectives are to select the elements of the 3x8 feedback matrix 
G so that the coefficient matrix of the closed-loop digital system, (p(T) - 6 (T)G, 
has the following form; 


<j)(T) - e(T)G 


( 6 - 8 ) 


where Aj Is 2 x 2 , A 2 ^3 s'"® 3x3 matrices, and the eigenvalues of Aj , A 2 and 
Aj are to be placed so that the bandwidth requirements described earlier are 
satisfied. 

Before the design can be carried out, one difficulty remains. It Is necessary 
to first calculate the matrices <J)(T) and 9(T), knowing A, B and the sampling period 
T. Since A Is 8 x 8 in the present case, the Laplace transform method becomes a 

I 

tedious task. The method chosen here is that of approximating (j>(T) by a truncated 


power series ; 1 .e. , 


<I>(T) = I AV/n! 
n=0 


(6-9) 


where N Is some positive Integer. 
The matrix 0 (t) is written 


e(T) = / (J>(T - T)dTB 

0 

T " A^(t -t)" 

= / I ^ 

0 n =*0 


(6-10) 


If the function 9 (t) Is uniformly convergent, the Integral and the summation In 
the last equation can be Interchanged. Also, truncating the Infinite series at 


n=N, we have 


e(T) = f J AV’^V(n+l)! B 
'^n=0 ^ 


( 6 - 11 ) 


For r = 0.02 sec, the results of (})(T) and 9(T) are given as follows: 



(j>(0.02)« 


9.998x10“' 

1.999x10"^ 

0 

-2.968x10*^ 

-1 .979x10*^ 

0 

2.968x10*^ 

1.979x10*^ 

-1 ,22l<xl0"^ 

9.998x10*' 

0 

-2.968x10*^ 

-2.968x10*^ 

0 

2.968x10*^ 

2.968x10*^ 

-4.1 60x1 0“^ 

-2.080x10*^ 

1 

1 .999x10*^ 

1 .999x10*^ 

0 

1 .011x10*^ 

5.059x10*'^ 

-6.240xl0“^ 

-4.160x10*^ 

0 

9.999x10*' 

1 .999x10*^ 

C 

1 .517x10*^ 

1 .011x10*^ 

-6.240xl0“^ 

-6,240x10*^ 

0 

-1.517x10*^ 

9.999x10*' 

0 

1.5)7x10*^ 

1.517x10*^ 

4.160x10“^ 

2.080x10*^ 

0 

1.011x10*^ 

5.059x10*'^ 

1 

1 .999x10*^ 

1 .999xlO*V 

6.240x10“^ 

4.160x10*^ 

0 

1.517x10*^ 

1.011x10*^ 

0 

9.999x10*' 

1 .999x10*^ 

6.240xl0“^ 

6.240x10*^ 

0 

1 .517x10*^ 

1 .517x10*^ 

0 

-1.537x10*^ 

9.999x10*' 


( 6 - 12 ) 


6 ( 0 . 02 ) = 


1.164x10*^ 

5.937x10*5 

-5.937x10*5 

1.164x10*^ 

5.937x10*5 

-5-937x10*5 

3.958x10*^ 

2.023x10*^ 

-2.023x10*^ 

5.937x10*^ 

3.035x10*5 

-3.035x10*5 

5.937x10*^ 

3.035x10*5 

-3.035x10*5 

-3.958x10*^ 

-2.023x10*^ 

2.050x10*^ 

-5.937x10*5 

-3.035x10*5 

3.075x10*5 

-5.937x10*5 

-3.035x10*5 

3.075x10*5 


(6-13) 


For the same bandwidth requirements described in Chapter 5, the following 
eigenvalues are selected in the s-piane: 

X component: (0.04 Hz bandwidth) 
p^, P2 “ -0.178 + j0.178 

(f), component: (10 Hz bandwidth) 


1 16 


P/^, P5 ■ '*^^•^36 + J^^.436 

(^>2 component; (I Hz bandwidth) 

Pg ” -20 

Py* Ps ° -^*^^36 + j4.A^36 

Using the transformation z = these eigenvalues are transformed Into 

the z-plane. The corresponding z-plane eigenvalues are: 

x; z,, Z2 » 0.996^^57 + JO. 0035^17 

<J)j : Zj = 0,0183156^} 

z^j, Zg = 0.25929^5 + JO. 31919^8 
<j>2: Zg = 0.67032 

Zy, Zg = 0,911366 + jO. 096987 

The condition of decoupling of the closed-loop digital system given In Eq. 
(6-8) results In kZ constraint equations. This Is due to the fact that Eq. (6-8) 
represents 64 scalar equations, but there are 42 zero elements In the matrix. 

The pole placements of the eighth-order system would produce 8 more constraint 
equations, so there are only 24 unknowns In the feedback matrix G. Therefore, In 
general, there would not be a set of solutions with more equations than unknowns. 
However, a closer look at the elements of iJ)(T) and 9(T) reveals that some of the 
elements are extremely small so that a good approximation can be obtained by 
assuming that these small elements are zeros. In other words, the state transition 
matrix <|>(0.02) given In Eq. (6-12) can be approximated by the following matrix. 
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9.998x10'’ 

1.999x10"^ 

1 ° 

0 

0 

1 

|0 

0 

0 

-1 .22i|xl0"^ 

9.998x10"’ 

’ 0 
f 

-2,968x10“^ 

-2.968x10'^ 

1 

|0 

2.968x10’^ 

2.968x10'^ 

0 

0 

L. 
! 1 

1.999x10'^ 

1 .999x10"^ 

r • 
10 

0 

0 

0 

0 

0 

9.999x10'’ 

1.999xlo'^ 

1 

1 0 

0 

0 

-6.2i|0xl0“5 

-i». 160x10“^ 

'1 0 
1 

-l.5l7xlo'^ 

9.999x10'’ 

1 

|0 

1 .517x10"^ 

1.517x10'^ 

0 

0 

r 

1 

0 

0 

T 
1 ' 

1 .999x1 0'^ 

1.999x10'^ 

0 

0 

1 

1 0 

0 

0 

1 

1 0 

9.999x10'’ 

1 .999x10'^ 

6.240x10'^ 

6.240x10“^ 

1 

10 

1.517x10'^ 

1.517x10"^ 

1 

1 0 

-1.537x10'^ 

9.999x10'’ 


(6-H) 

Similarly, we can approximate 6(0.02) by the following matrix: 



0 

0 

0 


1 .164x10'^ 

5.937x10"^ 

-5.937xlo'^ 


0 

0 

0 

0(O.O2)is 

0 

5.937x10“^ 

0 

3.035x10"^ 

0 

-3.035x10'^ 


0 

0 

0 


0 

0 

0 


-5.937x10"^ 

-3.035x10'^ 

3.075x10“^ 


Then, the matrix 0(O.O2)G becomes 


OF POOR 


(6-15) 


I 



e(0.02)G s: 


1 18 


“ 0 

0 

1 0 
) 

0 

0 

1 0 
I 

0 

0 

(66)2, 

( 60)22 

1 (00)23 

(00)2/, 

(66)25 

1[60>26_ 

(00)27 

( 60)23 

0 

1 

0 1 0 

0 

0 

0 

0 

0 

0 

0 

1 ° 

0 

0 

1 ° 

0 

0 

(66)5, 

(00)32 

1 

1 (00)53 

(60)5^ 

(00)55 

!( 6“>56 

T“ •“ - - 

(68)52 

(66)58 

0 

0 

1 0 
1 

0 

0 

^ 0 
1 

0 

0 

0 

0 

1 ° 

0 

0 

1 

1 0 

0 

0 

_(6G)g, 

(60)32 

1 <®<=>83 

( 60 ) 3 ^ 

( 60 ) 35 

1 

1 (66>86 

(60)37 

(66)88_ 


(6-16) 
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Based on Eqs. (6-1^), (6-15) and (6-16), the following 16 equations are 
obtained for decoupling the (|) - 0G matrix: 

(00)23 ^21^13“^ ®22^23 ®23^33 “ ^23 “ ° 

(BG)^/, = Qei^l't ®22^2i« ®23^3if “ ‘*’24 ’ “2.968x10 ^ 

(eG>25 - 62,9,5 + 622625 + 623635 “ ♦25 ■ -2-S68XI0-7 

(0G)2g = 021^16 ■*“ ®22^26 ®23^36 “ ‘*’26 “ ^ 

(88)2^ “ ^21^17 ^22^27 ^23^37 ™ ‘*’27 ~ 2.968x10 

~ ®21®18 ®22®28 ®23^38 " '*’28 “ 2.968x10 

(0G)gl ~ ^51^11 ^ ®52®2l ^ ^53^31 ~ ‘*’51 ~ "6.240x10 

(0G)j2 “ ®5i9i 2 ®52^22 ®53^32 “ ‘*’52 ~ "4.160x10 ^ 


(0G)j^ = 95|9 i6 + ®52926 ®53^36 " '*’56 “ ° 

(00)^^ = 05,g|y + 652927 ^ ®53^37 " %7 “ ^-517x10 


-5 


(0G)gg - 05]g|3 + 652928 ^53638 " ‘*’58 ■ ' * 517 x 10 

(0G)g^ ” ®8 i6i1 ® 82®21 ®83^32 “ '*’81 “ 6 . 240 x 10 


-7 



119 


(0G) 


“ 6.240x1 o”® 


( 00 ) 

( 00 ) 

(00) 


83 


84 


^'82 
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.517x10 


.517x10 


"5 


82 “ ®8)3i I ^ ®82^22 ®83®32 

®81^13 ®82^23 ■** ®83^33 

®81^14 ®82®24 ®83^34 “ ‘*’84 

'85 " ®8i9i5 ®82^25 *** ®83®35 “ ‘*’85 

The eight additional constraint equations for the solution of the 24 
elements of the feedback gain matrix G come from the pole placement requirements 
on the submatrices, Aj , A2 and A^. 

After the matrix <j) - 0G has been decoupled, It can be written as 


<J)(0.02) - 0{O.O2)G = 


where 


A, = 


‘*’11 

‘*’21 



0 

0 1 

_ _ 1 

0 1 

Ao 

0 

- - J- 



1 

0 1 


‘*’ 12 ' 


0 

‘*’22 _ 


_(eG)j 


22 


(6-18) 


‘*’11 

^2] - (9<5) 


21 


'*’12 

<^22 - 


22 


(6-19) 



<l>33 

‘*’34 

* 35 ' 


0 

0 

0 

''2 = 

^3 

‘*’44 

'•''15 

- 

0 

0 

0 


^<j>53 

‘*’54 

♦ 55 . 


(60)53 

CD 

(00)55 

r 

<(>33 


‘*’34 

■>35 

1 



‘*>43; 

<t>53 - (00)53 '*’54 " 


54 ‘*’55 


‘*’45 

( 00 ) 


55 


( 6 - 20 ) 
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‘*’66 

‘*’67 

‘*’68 


0 0 

A ^ 

A3 

‘*’76 

'*’77 

‘*’78 

- 

0 0 


_^86 

*e* 

CD 

'*’88 


(92)00 (02)37 

r 

"‘*’66 


‘t’67 


'*’68 


‘*’76 


'*‘77 


‘*’78 


‘*’86 " 

(ee)s6 

(*’07 - (00)37 

col 

CO 

0 

o> 

1 

CO 

CO 


0 

0 





( 6 - 21 ) 


The characteristic equation of the decoupled x dynamics Is 


z - (^1, - <|),2 

I zl “ A I ^ 

’ ■* '^>21 ■** ^®^^21 ^ “ '^'22 ^®^^22 

= + (( 92)22 " '{’11 " ' t ' 22 )^ *^ n '^22 ' ‘ f ’ ll ^®^^22 ■** ‘^’ 12(^®°^21 " *^21 

= (z + aj+Jbp(z + 3j-Jb^) (6-122) 

where 3| = -0.99641»57 and bj = -0.0035^17* 

Now solving for the unknowns in Eq, (6-22), we have 


( 92)21 = 021^11 ®22^21 ®23^31 

( 92)22 “ 021^12 ®22^22 '*' ®23^32 “ ‘^'l 1 ^1^12 "** ^^1 


(6-23) 
(6-2A) - 


For the <{)j dynamics, the characteristic equation of A 2 Is 


z - <l>33 “ ‘}>34 “ <*’35 

^ "<*>/{ 3 ^ " '*’44 " ‘*’45 

- (J)53 + (00)53 ■ + (00)5/, z “ <l>55 + (00)55 

I 

= + ((00)55 “ <l>55 ■ <l>z,4 " '*’33) ^■^‘*’44 '*’33) '*’45^^*'^ 54 
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**' 53 " '^skhs * ‘^33*^55 ‘**33‘^^t^i ’ " '^35*^53 * ‘^^^‘*’45-’ 

+ (<‘J>33‘I>4/, " ‘f'43‘J’3/,H®®)55 ^‘*’35‘*’43 " ‘<’ 33 ^ 5 ^ 

+ (<l>3/j<l>/,5 - <j>35'l>/,z,)(0G)53 " <^>33'<»/»i^‘^'55 + ‘{»33%i,'J>i,5 

“ ‘^43‘^’35‘^5'» " '^35'^A4%3) 


The desired characteristic equation Is written as 

(a ;• 32 + Jb2)(z + ^2 - Jb2)(z + C2) “ 0 


(6-26) 


where 


C2 » -0.01831564 
32 “ -0.2592945 

62 “ -0.3191948 


Equating Eqs. (6-25) and (6-26), we get 


(00)55 “ ®51^15 ■^ ^52^25 ®53^35 

(^*^^54 “ ®5i9)4 V24 •*• ®53^34 


<l>55 + + 1 + 232 "^2 


(6-27) 


•J, (I + 232 ^ I 


' *45(^35 - - ^35*w) " nF * '”*55 * 

+ 1 t 2aj + Cj) - + *^5*5^ - 1I.55 - + 2ajCj + aj + bj) (6-28) 


(00)55 “ ®51^13 ®52^23 ■*■ ®53^33 

(I + 232 ^ ®2 ^ ^2^^1 ■*■ '^2^ 


‘^‘35 ‘^34‘<’45 ■ '^35‘*’44 

For the ^2 dynamics, the characteristic equation of Is 
^ ■ *66 ■ *67 


(6-29) 


zl 


- <1, 


76 


<J> 


77 


<j)86 + ^ 


86 


^87 


87 


'88 


88 
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- z + ((0G)88 ■ <^>88 " *^77 " ■** (“^'^77 + 88 * *^78^®*^^ 87 

+ <|>68(0G)8g - <{.37^73 + <j)g^<j,Qg + + ■}>77«73 )z 

(('i'66‘^’77 “ '^76*^67 ^^®®^ 88 ^hd'^76 ~ '<> 66 ‘*’ 78 ^ ^®®^87 

■*■ ^'<’67‘^78 " ‘•’68‘**77^ ^®®^86 “ ‘^66'*’77‘^88 ■*' h6'^87^78 ’*’ ^7Sh7'^88 
“ •<»76‘^68‘<’87 " ‘^67‘*’78‘^86 '<' 68 ‘<’ 77 ‘<’ 86 ) ( 6 - 30 ) 

The desired characteristic equation Is written as 

(z + a^ + Jbj)(z + “Jb 2)(2 + c^) ■ 0 ( 6 - 31 ) 

where » - 0.67032 

ag » -0.911366 
bg » -0.096987 

Equating Eqs. (6-3O) and ( 6 - 31 ), we get 
(®®^86 “ ®8i9i6 ®82926 * ®83936 

(1 + 2a^ + a^ + bj) (1 + c^) 

•^68 ‘*’67‘**78 " ‘I’68‘*’77 (8-32) 

( 9 G)g 7 = 0 g,gj 7 + 032927 + 833937 

<)) 38 (l + 233 + 83 + b 3 )(| -f C 3 ) I 

~'J’78^‘<’68 ‘*’ 67‘^78 " hd‘^77^ ^ ^*■^‘^77 ■*■ ^*^88 ■*“ ‘*’77 ■** ' + 2 a 3 +c 3 ) 

- ‘^’78*^87 " ‘^’88 ■ ‘^’77 ^3 '’3) ( 8 - 33 ) 

^®®^88 “ ®8)9|8 ®82928 ®83938 

= <I >88 + 'I '77 + I + 283 + C 3 ( 6 . 3 /,) 

Equations (6-23), (6-24), ( 6 - 27 ), (6-28), (6-29). (6-32). (6-33) and (6-34) 
form the additional eight constraint equations, together with the 16 equations 
In Eq. (6-17) are adequate for the solution of the 24 unknown feedback gains In G. 
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The solution of the Zk constraint equations gives the following results; 
g,l - 35 .^ 3251359590 I »985 

9 i, - 203 . 57 HI 8')7387285 i>aGE IS 

gjj - -9.3684811250209332x10' Ol?' ^ 

g^j - - 8 . 11397701621271032 x 10 ^ 
g ,5 - - 1 . 4453443216549823 x 10 ^ 

9l6“9l7"9l8"° 

g^, “ “ 71 . 36)744259885901 

®22 " “^ 09.94846384061196 
g^j » 2 . 2392641102247482 x 10 ^ 

» 1 . 9394105849068884 x 10 ^ 
g^g » 3 . 454675889387745 x 10 ^ 

926 ** 3 . 578272520216161 x 10 ^ 

9^7 = 1 . 3239693135355619 x 10 ^ 

928 “ 1 •2749737170875715x10'* 

931 “ 932 “ 0 

9gg » 4.015292755500438x10^ 

93 /, = 3 . 4776162570406019 x 10 ^ 

9 gg = 6 . 194684677849835 x 1 o'* 

936 “ 3 . 578272527392904 x 10 ^ 

937 » 1 . 323969309412944 x 10 ^ 

9 gg = 1 . 27497371 29324844 x 1 o'* 


